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We derive perturbatively the gap equations for a color-superconducting condensate with 
total spin J = in dense QCD. At zero temperature, we confirm the results of Son for 
the dependence of the condensate on the coupling constant, and compute the prefactor to 
(^ ' leading logarithmic accuracy. At nonzero temperature, we find that to leading order in 

^Nj , weak coupling, the temperature dependence of the condensate is identical to that in BCS- 

like theories. The condensates for total spin J — 1 are classified; to leading logarithmic 
accuracy these condensates are of the same order as those of spin J = 0. 
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Cooper's theorem |Q-^ implies that if there is an attractive interaction in a cold Fermi sea, the system is 
0^ ■ unstable with respect to formation of a particle-particle condensate. In QCD, single-gluon exchange between 

quarks of different color generates an attractive interaction in the color-antitriplet channel |j]. Thus, it 
appears unavoidable that color superconductivity occurs in dense quark matter which is sufficiently cold 
p-pO[. How a dense quark phase matches onto hadronic matter is difficult to address R,01- In particular, 
while a quark-quark condensate may form, such condensation competes with the tendency of a quark-quark 

^ ' pair to bind with a third quark, to form a color-singlet hadron. 

r^ , One way of understanding color superconductivity is to compute at very high densities, where by asymp- 

totic freedom, perturbation theory can be used. At nonzero temperature, but zero quark density, it is well 
known that perturbation theory is a particularly bad approximation |21| . If 5 is the coupling constant for 
QCD, the free energy is not an expansion in g^, but only in g, with a series which is well behaved only 

V^ I for g < I. In contrast, at zero temperature and nonzero quark density, the free energy is an expansion in 

g"^ ln(l/ff) and appears to be well behaved for much larger values of the coupling constant, up to values of 
g < 4 |g2[ . Similar conclusions can be reached by comparing the gluon "mass" , nig ~ g/i or ~ gT to the 
chemical potential, fi, or the temperature, T [ |l5| . Thus for cold, dense quark matter, perturbation theory 
might give us information which not even lattice QCD calculations can provide. 

Color superconductivity is rather different from ordinary superconductivity, as in the model of Bardeen, 
Cooper, and Schrieffer (BCS) R-H]. In BCS-like theories, superconductivity is determined by infrared 
divergences which arise in the scattering between two fermions close to the Fermi surface: the initial fermions, 
with momenta k and — k, scatter into a pair with momenta k' and — k'. Summing up bubble diagrams 
generates an instability which is only cured by a fermion-fermion condensate. If the fermions interact 
through a point-like four-fcrmion coupling, though, there is no correlation between the initial and outgoing 
momenta, k and k'. In the gap equation, this implies that the gap function is constant with respect to 
momentum, as long as the momenta are near the Fermi surface. 

In QCD, however, scattering through single-gluon exchange strongly correlates the direction of the in- 
and out-going quarks: there is a logarithmic divergence for forward-angle scattering, ^ J d6/0. This extra 
logarithm from forward-scattering implies that the gap is not an exponential in l/.g^, as in BCS-like theories, 
but only in 1/g |7P, [l^ , [l5| . As a consequence, the gap function is no longer constant as a function of 
momentum, even about the Fermi surface. 



The logarithmic divergence for forward-angle scattering arises because in cold, dense quark matter, static, 
magnetic interactions are not screened through a "magnetic mass" . This is very different from a system of 
hot quarks and gluons. Over large distances, a hot system is essentially three-dimensional; gluons in three 
dimensions have power-like infrared divergences, which screen static, magnetic fluctuations through a mag- 
netic mass ~ g^T . In contrast, loop corrections in cold, dense quark matter are essentially four-dimensional; 
infrared divergences are at worst logarithmic, so that any magnetic mass is at best ~ /i exp(— l/g~), which 
is much smaller than the scale for color superconductivity, ^ /i exp(— 1/5) WJq]. 

The dependence of the zero-temperature color-superconducting spin-zero condensate 0o on the QCD 
coupling constant g was first computed by Son S, who used a beautiful renormalization-group analysis to 
show that 

-^o = 2^^exp(^-^j . (1) 

Here g — g/(3v27r) arises naturally from the solution, assuming three colors. In Son's result, 2 60 is a pure 
number of order one. 

In this paper we derive the gap equations for a condensate with total spin J = and Nf — 2 massless 
quark flavors, at an arbitrary temperature T. Our derivation is perturbative, i.e., 5 <C 1. In this case 
there are three scales in cold, dense QCD, the chemical potential fi, the gluon mass nig ^ gfj,, and the 
color-superconducting condensate 0o ~ M exp(— l/g), and they are naturally ordered, fi ^ nig ^ (J>q. 

We solve the gap equations to "leading logarithmic accuracy" , by which we mean the following. In the 
gap equations, the leading terms are ~ ln(/Lt/^o)- These terms generate the exponential in 1/g in Eq. (|l|), 
therefore ln(/i/(/)o) is of order 1/g. There are also leading logarithmic terms ~ In^fi/nig), which are ~ ln(l/g) 
plus a constant. The ln{l/g) gives rise to the prefactor l/g^, while the constant contributes to bo. For Nf 
flavors of massless quarks wc find 
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6o = 256^^(^— j b'o. (2) 

There are other terms in the gap equations, which do not arise from ln(/i/rng). These terms are of order 
one, and thus of the same order as the constant term originating from the logarithm hi(fj./mg). Hence they 
contribute in the same way to 69. We do not compute these, so that in (0) there is an undetermined constant 

Our results were described previously IQ , and agree completely with an independent analysis by Schafer 
and Wilczek pSJ ; they also overlap with those of Hong et al. [E2| . The factor Nr in Eq. (0) originates 
from the Nj dependence of nig-, however, 6g will also depend on Nf JToj . 

Even though we do not determine Vq, it is very interesting that the numerical value of b o/b[^ i s large. This 
implies that, for chemical potentials of order ~ 1 GeV, the gap can be of order 100 MeV 0,ll^. Such large 
values of the gap are in accord with previous estimates obtained within Nambu-Jona-Lasinio models |q], 
and are much larger than original estimates by Bailin and Love, 0o ~ 1 MeV [Q. 

We then solve the gap equations at nonzero temperature. We find a surprising result: while the detailed 
form of the gap function is very different in QCD versus BCS-like theories, the temperature dependence of 
the condensate — the dimcnsionless ratio of the condensate at a temperature T to that at zero temperature, 
<j){T)/ (j)^, is identical to BCS-type theories. In particular, the ratio of the critical temperature, Tc/4>o — 0.567, 
is as in BCS [|l|,|[. Our result for (f){T)/(f)o is valid to leading order in weak couphng, even though we cannot 
compute the overall magnitude of (pQ, i.e., b'g. In our mean- field approximation the transition is of second 
order, but the transition can be driven first order by critical fluctuations near the would-be critical point M. 

We then classify the condensates with total spin J = 1. Following our classification of condensates with 
J — 0, which employed projectors for chirality, helicity, and energy, we show that there are two types of 
spin-one condensates, longitudinal and transverse. As an example, we solve the gap equations for Nf = 1 
and find 

^1 =2^/iexp (-— ) , bi = -^b[. (3) 

r V 2gy 60 

That is, not only is the parametric dependence of the spin-one condensates on g the same as for spin zero, as 
argued originally by Son Isj, but to leading logarithmic accuracy, even the constant in front is the same. We 



do not expect that b'^, the undetermined constant analogous to b'^, is the same. Surely b[ is smaller than 6g, 
and depends on whether the condensate is longitudinal or transverse. This is very different from BCS-like 
theories, where condensates with higher spin are typically exponentially suppressed. In QCD, higher-spin 
condensates are only suppressed by a pure number. 

Nevertheless, one dramatic implication of our results is that superconductivity in QCD may be very unlike 
one's intuition from non-relativistic systems. Instead of higher-spin gaps being much smaller, they may be 
relatively large, which is important for phenomenology. Consider, for instance, a quark star with u, d, and 
s quarks. In the limit of high densities, when the strange quark mass rris is negligible, the number of u, 
d, and s quarks are equal. Then charge neutrality is automatic, and the preferred color-superconducting 
condensates are of spin J = 0, with color-flavor locking Q. For realistic densities, however, the chemical 
potentials for u, d, and s will not be equal. The strange quark chemical potential differs from the up and 
down quark chemical potentials due to Ws ^ tti^, nid- The up and down quark chemical potentials differ 
on account of charge neutrality. These effects suppress the color-antitriplet J = condensates, because they 
are composed of quarks with different flavor. On the other hand, the color-antitriplet J — 1 condensates 
may form between quarks of the same flavor, and are not suppressed if the chemical potentials of the various 
flavors differ. 

This paper is organized as follows. In Section O, we derive the gap equations for a spin J — condensate 
of Nf — 2 massless flavors at an arbitrary temperature T. In Section |l| we solve these equations to 
leading logarithmic accuracy, first at T = 0, and then at nonzero T . In Section |^ we classify the possible 
spin J — 1 condensates, and solve the gap equations for Nf — 1. Section ^ concludes this work with a 
discussion of higher-order effects which contribute to &q and b'^. Our units are ?i = c = 1, the metric tensor 
is g^^'' = diag (-|-, — , — , — ). 4-vectors are denoted by capital letters, K = K^ = (A;°, k), and k = |k|. 

II. THE GAP EQUATIONS 

In general, a color-superconducting condensate ^J- r is a Nc x Nc matrix in fundamental color space 
{i,j = l,..., Nc), a Nf X Nf matrix in flavor space {f,g — 1, . . . , Nf), and a 4 x 4 matrix in Dirac space |0]. 
As shown in |0], for Nf = 2 a color-antitriplet condensate is a flavor singlet, 

Kj.fg = ^f9 *y- = e/9 ^Vk '^t ■ (4) 

By a global color rotation, we can always choose ^^ to point in 3-direction in color space, ^^ = Sks $^. 
For Nf — 3, even if we assume that the dominant condensate is a color antitriplet, there is always a small 
admixture of a color sextet |lq|. The color-antitriplet condensate is a flavor antitriplet, and the color sextet 
a flavor sextet [[7|jl9|]. In this paper, for simplicity we consider only the case Nf = 2. 

The gap equation for a color-superconducting condensate of massless fermions was derived in Q [cf. Eq. 
(A35) of H]]. At a nonzero temperature T, this equation reads (we suppress fundamental color and flavor 
indices for the moment): 

^+iK)=g'^Y.^^a^ll{K~Q)G^{Q)<^+{Q)G+{Q)V^ . (5) 

Q 

Here, T/VJ2q = '^J2nl^'^^/i'^'^)^ i^ ^^^ infinite-volume limit {n labels the Matsubara frequencies ojn ^ 
(2n -|- l)7rT = iqo), and a summation over Lorentz indices fijV as well as adjoint color indices a,b — 
1, . . . , Nc — 1 is implied; A^c = 3 is the number of colors. A°^ is the gluon propagator, 

G±(Q)^(7.Q±M7o)-^ (6) 

is the propagator for free, massless particles (upper sign) or charge-conjugate particles (lower sign), 

is the propagator for quasiparticles (upper sign) or charge-conjugate quasiparticles (lower sign), and 

I]± = $=FGT$± (8) 



their self energy arising from the interaction with the condensate. The charge-conjugate condensate is 

$-=7o($+)So, (9) 

and the vertices are 

r^^T,7^ , f^ ^ c (r^f c-1 ^ -7''Tj , (10) 

where 7'' are the Dirac matrices and Ta the Gell-Mann matrices, C = — C^ = —C'^ = —C~^ = 17^70 is the 
charge conjugation matrix. 

In the following, we analyze the flavor, color, and Dirac structure of the gap equation (^. For Nj = 2, 
the color-flavor structure of the condensate (0) does not mix color and flavor indices, and thus the analysis 
of flavor and color can be done separately. This is different for Nf = 3, where color rotations are locked to 
flavor rotations. 

A. Flavor structure 

We first discuss the flavor structure of Eq. (|^). Fundamental color indices will be suppressed for the 
moment. The free propagator (|6|) is diagonal in flavor, 

Go%g(Q) = ^/s G±(Q) ^6f,ij-Q± ^,Jor' ■ (11) 

So is the self energy, 

^/9 = 70 (<i>+)^,, 70 G-^, $+ ^ e^f S^, e.g 70 ($+)So G^ $+ = Sj, 70 ($+)So G^ cO+ ^ 5fg S+ , (12) 

and thus Gt — Sfg G+. Therefore, both the left- and right-hand sides of (||) are simply proportional to e/gi 
the flavor structure of the gap equation is trivial in QCD with Nf ~ 2 flavors, and will thus not be explicitly 
denoted in the following. 

B. Color structure 

The free propagator (|g) is diagonal in the color indices for the fundamental representation, 

Go%(Q) ^ % GtiQ) ^ S,, (7 • Q ± A^7o)"' . (13) 

The self energy T,f- is also diagonal, but not all diagonal elements are equal: 

S+ = 70 (*+)L 70 G,,i $+ = ek^3 Ski ei,3 7o (*+)So Go" $+ ^ (<5., - ^,3 ^,3) S+ (14) 

— the self energy for quarks with color 3 vanishes ! This is easy to understand. Let us first note that, 
according to Eq. (A17) of |^, the condensate <&+ is actually proportional to (V'cV'), while $~ '^ (V'V'c)- 
(ipc = CtJ)'^ is the charge-conjugate fermion field.) Thus, according to Eq. (g]), the self energy S+ arises 
from the following process: a quark with, let's say, color 1 annihilates with a corresponding antiquark in 
$"*", creating a charge-conjugate quark with color 2. This quark is propagated with the charge-conjugate 
propagator G'^ , and annihilates with a charge-conjugate antiquark of color 2 in $^, whereby a quark with 
color 1 is emitted. As only quarks with colors 1 and 2 condense, it is not possible to annihilate and create 
quarks with color 3 in this process; thus the latter do not attain a self energy. 
One can now compute the color structure of the quasiparticle propagator, 

G+- - {% [G+] "' - (% - <5,3 'Jjs) E+}"' = (5„ - 5a <5j3) G+ + ^,3 ^,3 G+ , (15) 

whereG+ = {[G([]"^-S+} . 



When inserted into the gap equation (J^), the terms ^ St^Sj^ in (|l5| ) vanish, as ^^ ^ eij^^. The gap 
equation becomes 

e,,3 $+(/0 = - {n n, T^^ Ty) 9% E ^^ ^^' (^^ - Q) Go"(Q) *+(Q) G+(Q) 7"^ . (16) 

Q 

From the exphcit form of the Gell-Mann matrices we now infer that only gluons with adjoint colors 1, 2, 3, 
and 8 participate in the gap equation. 

A color-superconducting condensate $^ ~ (5^3 breaks SU{i)c to SU{2)c] gluons 1, 2, and 3 then correspond 
to the generators of the unbroken subgroup, and thus remain massless, while the 8th gluon attains a mass 
through the Anderson-Higgs effect. Denoting A^^ ~ A^^ = A|^^ = A^i, and A^^ = A^^, we obtain 



^-^w-I^^^E 



4^ V ' ' ' 



1 ~. 
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^^,,{K-Q)--K^,{K~Q) 



Go-(Q)$+(Q)G+(Q)7% (17) 



where a common factor e^js has been dropped from both sides of the gap equation. 

In a complete treatment of the gap equations the effect of the condensate on the gluon propagator has to 
be included [Q. In contrast, we use the gluon propagator in the "hard dense loop" (HDL) approximation 
p3|-p6[. The HDL propagator introduces a gluon mass rUg ~ 5/i. Since in perturbation theory the scale of 
the condensate 0o is much smaller than nig, it is reasonable to expect that, to first approximation, we can 
neglect the effects of the condensate on the gluon propagator. A more detailed explanation is given below. 
We consequently assume A^i^ = A^j^, and obtain 

<^+{K)^\g^yY.^^^^^-^^-Q)G^{Q)'^+{Q)G+{Q)^^ . (18) 

■^ ^ Q 



C. Dirac structure 

In H we have shown that the gap matrix $+ for a condensate with total spin J = has the form 

J = : <i>+{K) = E E E ^ls{K)VU^) . (19) 

h=r,e s=± e=± 

The 4x4 matrices 7'^j,(k) are defined as 

VUk)^VhVsik)A'^ik) , (20) 

where Vh are projectors for chirality, h — r,(,, 

Vr^'-V^, V,^'-^, (21) 



while 



are projectors for helicity, and 



^ ,, ^ 1 + s 75 70 7 • k , , 

>^s(k) ^ , s = ± , (22) 



^e(k) = i±£(Ai7o^Ji+i^k7o) , e = ± , (23) 



are projectors for energy, /3k = k/E^, ak = m/E^, Ey^ = y/k^ + m?. While the Vf^gik) are composed of 
projectors, they are not projectors themselves [cf. Eq. (B29) in Q], and thus were termed "quasiprojectors" 

in i- 

In the ultrarelativistic limit, to = 0, the energy projectors simplify to 



m^O: A^],) = l±£JpL^. (24) 



This has three major consequences. First, the quasiprojectors (20) become true projectors [cf. Eq. (B29) of 
y|]. Second, V~_^ = V^__ = 'P^_^ = 'P^_ = 0, expressing the fact that right-handed, positive- hehcity particles 
cannot have negative energy, etc. [cf. Eq. (B30) of 0]. Third, either the chiraUty, or the hehcity, or the 
energy projector in ( po| ) becomes superfluous [cf. Eq. (B31) of [Q|]. In the following, we omit the helicity 
projector and use just chirality and energy projectors, 

m - : Vi, (k) — > Vt (k) ^ Vh A" (k) . (25) 

Then, (^ simplifies to [cf. Eq. (8) of g|] 

J = 0, m - : <P+iK) = J2 H '^'^ W ^'"^ W . (26) 

h=r,£ e=± 

The quasiparticle propagator assumes the form [cf. Eq. (15) of [g|] 

G^(Q) = E E , 7^1,.2 [GoiQT' , (27) 



where ig{(f>) = \/{q — efiY + |0|^- From Eq. (26) of [|9[| we then infer that the gap equation (|l8| ) can be 
written in the form 

Q ;i=r,f e=± 90 ~ LSv^'i''J 

where —h = £, if h — r, and —h — r, if h — £. 

With the help of the projectors 7^^(k) one can derive gap equations for the individual gap functions (pf^, 

K^'iQ) .Tr[V^^{k)j^Vl^{<i)Y]] . (29) 



^—G ( J,— e\ 



To obtain this result we have used Vh 1^ — 7^ V-h and Vh V-h — 0, so that the gap equations for right- and 
left-handed condensates decouple. This is a consequence of the U{1)a symmetry of the QCD Lagrangian, 
which is expected to be effectively restored at asymptotically high densities ||7|,|l4[|. 
In Coulomb gauge, the HDL propagator is p^ , p4[ | 

Aoo(P) = AKP) + ^c^ , Ao.(P) = Cc ^ , ^^AP) = (-5., - PrP,) MP) + Cc ^ . (30) 

We set the gauge parameter ^c = 0. We shall show later that our results are manifestly gauge-invariant 
to leading logarithmic accuracy. The longitudinal and transverse propagators Aj.t are defined in Eqs. ( p3| ) 
below. 

The traces in ( p9|) are readily evaluated. We need the terms 

Tr [Vh A^(k) 70 AT^(q) 70] = ^^^ , (31a) 

J2 Tr [Ph A'^(k) 7. A^'=(q) 7.] = - ^^^ , (31b) 

i 

Tr[P,A^(k)7-pA^^(q)7-p] = _ i^iilS (^^ , (31c) 



where p = k — q. Obviously, the final result is independent of the chirality projector. We therefore conclude 
that the gap equations (p9) for right- and left-handed gap functions are identical. This means that right- 
and left-handed gaps are equal up to a complex phase factor exp(j6'). Condensation fixes the value of and 
breaks U{1)a spontaneously. As discussed in [0,H this leads to spontaneous breaking of parity. 
The gap equations for either right- or left-handed gap functions read 
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9o- h'(<^/i'). 






+ ^*^^ ^)|^ 2 ' 2 (k-q)2 


J . 



(32) 



The gap equations involve singularities from both the quark and gluon propagators. The poles of 
^/ \la ~ [^qi'l^h)] \ Si'^^ ^ residue '-^ l/e^((/)|). Remember, though, that the quasiparticle energy e+ is 
very small near the Fermi surface, e+ ~ \4'h\- The quasi-antiparticle energy e^, however, is always larger 
than ^. Therefore, in weak couphng, the dominant terms arise from the quasiparticle poles qo = ±e+(0^) 
[g[, and the contribution from quasi-antiparticle poles can be neglected. Consequently, we do not need to 
solve Eq. ( p2[ ) for the quasi-antiparticle gaps, (j)'^, in order to determine the solution for the quasiparticle 
gaps, 0^. In the following, we drop the subscript h and superscript + to simplify the notation, and denote 



D. Spectral representations 

To perform the Matsubara sum over quark energies qo = —i{2n+ l)7rr, we introduce spectral representa- 
tions. For the gluon propagators P3,e3, 



Ai{P) = -^+ f dTeP°^Ai{T,p) , MP)= f dTeP°^At{T,p 
P Jo JO 

/•oc 

Ai,t{T,p)= / dujpu{iJ,p){[l + nB{uj/T)]e~^' + nB{io/T)e^^} , 
Jo 



(33a) 
(33b) 



where nsix) = l/(e^ — 1) is the Bose-Einstein distribution function. The term —l/p^ in the longitudinal 
propagator cancels the contribution of A;(P) at po — > cx) |24|. The spectral densities are given by P3,M 



p/,t (w, p) = pf°t°{^, p)5[uj- uji^tiv)] + pUH^, p) &{p - ^) 

pole/ N Cj(w2_ 2-) 

Pi (w, p) = 

pr(^,p)- 



p2 (p2 + 3^2 _ ^2 


) ' 












2Af2 oj j 


p^ + 3m^(l-gln 


P + UJ 

p— CO 
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'-M) 


TT p y 


LoiiU^-p^) 




imloj-^-{uj-^-p^Y ' 


A-P LO p2 1 


\ 2 3 2 ( W' 

n 1 m 


f — In 
2p 


P + LU 
p — UJ 


)] 


TT p p2_ 


-C.2 


[P 1 2 ""« \p^ 


~UJ^ 



M^- 



(34a) 
(34b) 

(34c) 
(34d) 

(34e) 



The basic parameter of the HDL propagators is the gluon mass 



? 9 

We also found it convenient to introduce 



iVc 



Nf\9'T^ 
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M2= ^m2~2.36m2 
4 y y 



The functions CLi;.f(p) are the solutions of the equations 

p^ + 3ml 1 
3 



2p 



1(1 i\ " 



2 I ^t / 2 2\ 1 

— {p - w J In 



2p 



W; f LOi +p 

In 



= 
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(35) 
(36) 

(37a) 
(37b) 



and define the dispersion relations for longitudinal and transverse gluons, respectively. They satisfy w; t(p) > 

We argue later that the phase-space region which dominates the gap integrals is the nearly static, small- 
momentum limit. The gluon energies are on the order of the gap, w ~ 0, while the gluon momenta p are 
much larger than w, but much smaller than mg. In this limit, uj <ti p <ti rUg, 

^ut 2M2 uj 1 ^^^ AP Lop 
Pi {^, P) , o , o ?^? ' Pt i^' P) Z ,..n .2 • (38) 



TT p [p"^ + 3mg) 
We also introduce a spectral representation for the quantity 



TT p6 + (M2tj)' 



s(g) 



il-[^M)Y 



/ dTe9«^S(r,q) 
Jo 



(39) 



Neglecting the singularities of (j){Q) in the complex qo plane, and assuming 4>{qo,(l) to be an even function 
of go , we find 



where 



/"OC 

S(r,q)= / dojp{oj,q){[l-np{oj/T)]e-'^^-nF{u;/T)e'^^} , 
Jo 

<^(w,q) 



Pi^,(i) 



2lo 



■5{uj - eg) 



<^qW > 



(40) 



(41) 



and npix) = l/(e^ -|- 1) is the Fermi-Dirac distribution function. 

By neglecting the singularities of (j){Q), the only contribution to the spectral representation of S((5) is 

from the poles of 1/ -^ (/q — [eg(0)] >, which generates the delta function in the spectral density mW). This 

forces the energy in the gap function (f>{qo, q) to lie on the quasiparticle mass shell, qo = eg(0). 
The Matsubara sums over go can now be computed as (p = k — q) 

T5:AKX-g)S(Q) = -^|^|-itanh(|,) A + ^°°d..p,(a;,p) 
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1 



1 



2 \2T / \ko + Lo + eg Kq — uj — eg Kq — uj + eg ko + lu — eg 



1 



1 



1 



1 



1 



- coth (—) , , 

2 V 2T / yko + LU + eg ko ~ lu — eg ko — lu + eq ko + lu ~ eg 



(42a) 



Tj2^tiK^Q)EiQ) = -^^^ rdLuptiLu,p) 

90 ^^« -^0 



1 



1 



-t- 



1 



2 \2T J yko + LU + eg ko — lu — eg ko — lu + eg ko + lu — eg 



- cothf— ") ( 

2 \2TJ \ko + Lu 



1 



1 



1 



LU + eg ko - LU - eg ko — lu + eg ko + lu - eg 
where we have made use of ( |4l] ) and e*^"/^ = —1, since fco = —i{2n + l)7rr. 



(42b) 



E. Analytic continuation 



In field theories, the only physical quantities are those on the mass shell, such as S-niatrix elements. In our 
case, this implies that we need the gap function not for Euclidean values of the energy, ko = —i{2n + 1)ttT, 
but on the quasiparticle mass shell, for real values of ko = ek- This is achieved by analytically continuing 
^0 ^ £/£ + *7? Pi- The analytic continuation of the gap equation ( p2[ ) introduces a term ii] in the energy 
denominators in (E^). As l/{x + ir]) = P{l/x) — iTT6{x) (where P stands for the principal value), this generates 
an imaginary part for the function ^{efc,k). As shown below, this imagina ry p art is down by g relative to 
the real part. This justifies neglecting the singularities of (j){Q) in deriving ([4^). In the following, implicitly 
we take the principal value of all energy denominators. Finally, since the gap function 0(A:o,k) does not 
depend on the orientation of the 3- vector k, but only on its magnitude, fc, cf. (32), the on-shell gap function 
is a function of k only, 4>{ek, k) = (t>{£k, k) = (t>k- 



F. Contribution from longitudinal gluons 



The contribution from longitudinal gluons is give n by Eq. (42a). In general, this equation can only be 
evaluated numerically, using the full expressions (pj) for the spectral densities. We can proceed analytically 
by noting that, due to the factor l/e^, the integral over momenta q is dominated by the region close to the 
Fermi surface. In the following, we shall therefore assume that the momenta k and q are close to the Fermi 
surface, such that e^ '^-^ e^ ^ mg. Now rewrite 



1 



1 



1 



1 



UJ + eq + ek uj + eq-ek uj - {eq + ek) lo - {eq - ek) 



1 



Cg +efc 



e<7 -Cfc 



e<7 + efc 



Eg - £/£ 



UJ + eq + ek uj + eq~ek lo - {eq + ek) uj - {eq - ek) 
and make use of the sum rule p3| , p4| 

Pi{uJ,p) 1 . ,„ . . ,„ > 1 



dcj ■ 



pz 



Ai{0,p) , Ai{0,p) 



p^ + 3to2 



(43) 



(44) 



to obtain 



TY^A,{K-Q)E{Q) = -'t^^[li.nl.(^\ 



90 



2eq 12 



\2TJ 



p^ + Sm^ 



Ji(k,q) 



ICi{k,q) 



(45) 



where 



Ji(k,q)= / dcj 



Pl{^,P) 



and 



£q +f-k 



€q + ek 



Cg — fife 



eg — Cfc 



(e^ + efc) UJ + eq + ek uj - {eq - ek) uj + eq ~ ek 



/C,(k,q)=/ duj pi{uj, p)-coth(^^j 



(46) 



UJ + eq + ek UJ + eq- ek uj - {eq + ek) uj - {eq - ek) 



(47) 



The integral jTj consists of two parts, one from the pole term in the spectral density, and one from the cut 
term. The pole term is 



:^r^k,q) = i. '^'-^' 



p2 p2 ^ 3j^2 _ ^2 



Eg + fife 



Eg + £k 



Cg - Cfe 



Eg - e/c 



uJi-{eq + ek) u>i + eq + ek uji - {eq - ek) uji + eq - ek 



(48) 



As oji > nig ^ e/c ^ Cg, we may expand the energy denominators around u;/, 



T-polc 



J'r"(k,q)^- 



1 ujf-p' A + el 



p^ p^ + Zvfvi — ujf ujf 



(49) 



This contribution is quadratic in Cq/toi ~ efe/w/ <C 1 and thus neghgible. 



The cut term is estimated using the approximate form of the spectral density (38), 



jr\KcL) 



2M^ 
■np 

2M^ 
■np 



1 



im?^' 



duj 



Cfe 



Cfc 



Cfe 



eg — Cfc 



Sm^)^ 



(eg + ek) hi 



uj- {eq + ek) 



Cfc 



uj- {eq- ek) 



Sq — ek 



p 



Eg + ek+p 



+ {eq - ek) In 



Cfe -P 



Cq - fife + p 



(50) 



This integral is proportional to factors eqztek, which parametrically cancel the prefactor l/e^ in (Hq). (This 
is true except for p — > 0, where these factors cancel when expanding the logarithms. However, the Jacobian 
of the angular integration in the gap equation, J d cos 9 ~ J dpp, suppresses this contribution.) It is thus 
neghgible compared to the term '^ 2/(p^ + Sm^) in (^). 

We now turn to evaluate /C; (k, q) . After expanding the energy denominators in (B7h around uji > nig 3> 
£k '^ iq, the pole part reads 



^polc 



(k,q)=^-^ 



9 9 

ujf - p-" 



p2 p2 



3?7lg 



UJt 



k-Hw)^'i 



(51) 



The range of temperatures of interest are limited by Tc^ the critical temperature for the onset of color 
superconductivity. As Tc ~ 0o <C rUg < uji, coth((jj;/2r) ~ 1. The factor eq cancels the prefactor in (Ha), 
such that ICf° '^ is negligible compared to the term ^ ^/{p^ + 3m?) in ([45|). 

To estimate the cut term, we use the approximate expression (pS), and employ a relation similar to (H3), 

111 1 



(52) 





u + eq + ek u + eq- ek lo - (eq + ek) uj ~ (eq - 


^ek) 






1 


eg + Efe 1 eq-efc ^ ^q + ek ^ 


eq-ek 






cj + e, + Efe u + eq- ek w - (e, + ek) 


w- (eg-efc)_ 


! 


to obtain 


^-HKc,)^~'''\,' ,^, Td^cothr-) 

np {p^ + iml)^ Ja 2 V2T/ 






X 


' eq + ek eq- ek eq + ek 


1 eg - efe 




u + eq + ek u + eq-ek u- [eq + ek) 


w- (eg-eA;)_ 



(53) 

q - tfcjj 

For cu of order T or larger, this contribution is small, since then coth(Li;/2T) ~ 1, and thus all terms are 
proportional to factors egicfe. For w <C T, however, the large (classical) occupation number density of gluons 
can enhance this contribution. In the classical limit, one approximates coth(a;/2T) ~ 2T/uj. Simultaneously, 
as the range of validity of this approximation is uj <^ T, the upper limit of the integral should now be 
replaced by p* = min{T,p). Reverting the step ([5^), we obtain 



jcr 



(k,q) 



2APT 



np (p2 + 3to|)2 



In 



p* + e„ + efe p* + eq 



efc 



p* 



(eq + ek) P* - (eg -efc) 



(54) 



For T ^ 0, we may expand the logarithm to show that this contribution is quadratically small in T. On the 
other hand, for T -^ Tc, close to the Fermi surface e^ ^ efc ~ ^ 0, and the logarithm can be expanded for 
p* ^ eq ± ek- The result is proportional to e^ ± efc, which cancels the prefactor l/e^, and thus suppresses 
this contribution. We shall therefore neglect it in the following. 
The final result for the longitudinal contribution is thus 



TY,MK-Q)^iQ) 



go 



(filS) i tanh (^) -^ 
2 \2tJ p2 



2e„ 



3m^ 



(55) 



This result could also have been obtained by simply taking the static limit of the longitudinal gluon propa- 
gator, Ai(0, p), on the left-hand side in (paj), and performing the Matsubara sum directly. We conclude that 
the contribution of static electric gluons dominates over that of non-static electric gluons. We note that, 
while the individual terms ( [l9[ ) and ( pl|) exhibit an apparent infrared divergent prefactor l/p^, the sum of 
all terms is infrared finite, as can be shown by computing the contribution of electric gluons without using 
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G. Contribution from transverse gluons 

The contribution from transverse gluons is written in a form similar to ( |4q ) 



T^A,(i^-Q)S(g) 



2e„ 



^tanh(g;) Jt(k,q)+/Ct(k,q) 



(56) 



where 



Jt(k,q) = / dujpt{uj,p) 
Jo 



1 



w + Eg + efc Lj + Eg - Efc w - (eg + e/c) w - (eg - efc) 



(57) 



and 



/■oo ^ 

/Ct (k,q) = / (iwpi(w,p) - coth 
Jo ^ 



\2TJ 



^ + e-q + <^k ^ + eg-efc w - (e, + efc) w - (e, - efc) 



(58) 



The function J7t consists of a pole and a cut term. In the first, we expand the energy denominators around 
^t > "nig 3> eg ^ efc, to obtain to leading order in eg ± efc: 



7-polc 



J'r"(k,q)c.4 



9 9 



3mW - (Cj2-p2)2 



(59) 



For the cut term, we again employ (pq) to obtain 



J'r'(k,q) 



APp 



duj ■ 



p6 + (M2c^)^ 



1 



w + eg + efc cj + eg - efc w - (eg + efc) cj - (eg - efc) 



(60) 



The to integral can be performed analytically. Denote a = p /M . Then 



JT*(k,q) 



•M2 



duj ■ 



LO 



LO'^ + a 
a 



1 



1 



1 



1 



cj + eg + efc cj + eg - efc cj - (eg + efc) cj - (eg - efc) 



a^ + (eg + efc)2 a^ + (e, - efc)^ 



arctan I — 



£k 



a^ + (eg + efc)2 



In 



p + eg + efc 



P - (eg + efc) 



eg — efc 



a^ + (eg - efe)^ 



In 



P + eg - efc 



P - (eg - efc) 



The terms proportional to eg ± efc are of higher order and can be neglected, such that 



jr'iK^) 



p^ + M\eq + ekf p^ + M^Cq-ekY 



2 f AP 

arctan 



(61) 



(62) 



The arctan cuts J7t off for momenta p larger than ~ M . To make further progress, in the following we replace 
the more complicated arctan with a simple 9 function cutoff, 2 arctan(M2/p2)/7r — ^ 9{AI — p). However, to 
retain consistency with the sum rule p3|,p4| 



duj 



P«(^,P) _ 1 



p^ 



(63) 



we also have to modify the result for the pole term (|59|). To enact the modification, note that j7t/2 is 
identical to the left-hand side of the sum rule (|6^), when we set eg = Cfc = in (|57|). Thus, the most simple 
choice is 



J'r''(k,q)^-0(p-Af), 

pZ 



(64) 
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since then, for eq = £k = 0, 



Ji(k,q)=J'r (k,q)+J'r'(k,q) 



p^ 



(65) 



The function JCt also consists of a pole and a cut term. In the first, we expand the energy denominators 
around ujt > rUg ^ eq r^ ek, 



-pole 



/Cr^(k,q)^- 



^t-p' i ^oth (^] i^ 



(66) 



As the factor e, cancels the prefactor l/e, in (p6D, this contribution is of higher order and thus negligible. 
The cut term is estimated with (p8) to be (as before, a = p^/A'P) 



^r(k,q) 



ttA'P 



duj 



CO 



;2 + a4 ^°*^ {^) 
1 



Cfc 



Efc w-(eq + efc) 



6.) 



(67) 



For w of order T or larger, coth(Li;/2T) ~ 1, similar steps that led to (|6l|) now show that the final result is 
proportional to eq ± ek, and thus of higher order. For w <C T, however. 



^r(k,q) 



pT r .. 1 




1 


1 1 


1 




ttM- 
AP 

■K 


'■ J„ ™C^2+a2 

r pT 

_p6 + M4(e, + efc)^ 


-In 


eg + ffc tj + 

P* + e, + e/c 

P* - (eg + efc) 


eq- ek uj- {eq + ek) 

pT 
'^ f' + A-neq-ekY "" 


uj-[eq- 

P*+eq- 

P* - (eg - 


-e.)J 
- efc 
-efc)J 



(68) 



where as before p* = min(T,p), and terms proportional to eq ± ek have been neglected in the last line. For 
small T, an expansion of the logarithms shows that this contribution is quadratic in T. On the other hand, 
for T -^ Tc, an expansion of the logarithms for p* ^ e^ ± ek shows that it is proportional to eqzLek', it thus 
parametrically cancels the prefactor 1/eq. We shall thus neglect this contribution in the following. 
Our final result for the transverse contribution is therefore 



T^At(A'-Q)S(Q) 



go 



^'«' "^^ I tanh (4^) \\e{p- M) + e{M - p) 



2e„ 



\2TJ 



+ 



p^ + M\eq + ekf p^ + M\eq-ek) 



(69) 



H. The gap equation 



Collecting the results from the previous subsections, and replacing the angular integration J d cos 9 by an 
integral over p = |k — q|, the gap equation ([3^) reads 



0fc 
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247r2/c 



M+>5 ^ 
dq - 

.-S e, 



tanh I —2- I (bg 



fe+g 



dp p 



k — q\ 



{k + qf 



3m2 



Akq 



pz. 



■ e{p - M) + e{M ~ p) 



p 



p^ + MHeq + ekY p^ + MHeq-ekY 



Akq 



{e-q^f 
4fcgp2 



(70) 



Here, we restricted the q integration to a region fi — d<q<fj, + S, S^fj, around the Fermi surface, since 
we assumed that eq ^ ek ^ rUg ^ gfi <ti fi- Therefore, we may set k ~ q ^ ^. Also, to leading order we may 
neglect p^ terms with respect to {k + g)2 terms. Then, the gap equation (|70|) simplifies to 



247r2 



M+<5 



fi — d 



^ t.nh (I^\ 



\2TJ 



In 



/4Ai2 



\3m2 



In 



AP 



I- 



AP 



(71) 
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The first term in brackets, ~ ln(4^^/3mg), arises from static electric gluons [proportional to 2/(p^ + Sm^) 
in Eq. <JI(\}], the second, ^ ln(4/i^/Af^), from non-static magnetic gluons [proportional to 9{p — M) in Eq. 
([70|)], and the last term, ~ ln(M^/|e^ — e^|), from soft. Landau-damped magnetic gluons [proportional to 
9(AI — p) in Eq. (uOh]. These terms can be combined to give 



g2 rs 



187r2 



i^'-(lf)^"(pr^)*.. 



where we exploited the symmetry of the integrand in (uW around the Fermi surface to restrict the q — /.i 
integration to the interval fJ- < q < i-J- + S, i.e., we only integrate over momenta q from the Fermi surface up 
to ijl + 6. Moreover, we have defined 

/ r, \ 5/2 

The temperature dependence of nig, cf. Eq. (pa), has been neglected, as T ~ 0o ^ "^g ~ 5M- 

The result ( fn^) for b agrees with the analysis of Schafer and Wilczek |l^. It does not agree with that 
of Hong et al. [|12|, because they neglected the contribution of electric gluons. Note that our gap equation 



( [72| ) differs from that of Son [g[ and of Schafer and Wilczek |13| in that we integrate over momenta q, while 
they integrate over Euclidean energies go- The energy dependence of our gap function is always given by the 
quasiparticle mass shell (t>{K) = 0(efc,k). One advantage of our approach is that the extension to nonzero 
temperature is immediate. 

From Eq. (|70|), the gluon momenta which dominate the contribution of nearly static, transverse gluons, 
are p^ ~ M'^{eq ± efe)^, or for energies e,, efc close to the Fermi surface, p ~ (m^^)^/'^. While these momenta 
are small relative to the gluon mass rrig, they are much larger than the condensate (j). This is why we can 
use an HDL propagator in the gap equations, neglecting the effects of the condensate. 

The approximations which lead to (u2|) are only valid to leading logarithmic accuracy, in that we neglect 
terms of order 1 relative to hi{p/mg) ~ ln(l/g). Possible contributions of order 1, which contribute to 5o in 
Eq. (0), are discussed in Section M. 

III. SOLVING THE GAP EQUATION 

The gap equatio n (172 ) is an integral equation for the function (/)/£ = 0(efe,k). In general, it can only be 
solved numerically [|l3|. In this section, we first discuss the parametric dependence of the solution to ([T^) on 
the QCD coupling constant. To this end, it is instructive to understand which terms in (uW determine the 
exponent and the prefactor in Eq. (|l|) . We then derive an approximate analytical solution by converting the 
integral equation (p2h into a differential equation. The solution is discussed in detail at zero and at nonzero 
temperature. 

A. Parametric dependence of the solution on the QCD coupling constant 

Let us introduce the variable 

(74) 



3\/27r ■ 
At T = 0, the general structure of the gap equation (^) is 



k -' 


f f'd{q-^i) 
2 io eg 


/ 
In 
\ 




) + 2 In 


f bo 



h ■ (75) 

The first term, ~ ln(^2/|e^ — e^|), arises from nearly static, transverse gluons. As compared to (|7l| ) we 
have factored out a term ~ h\{M'^ /fj?). The second term, ^ ln[&o/(&off^)]i combines this term with the 
contribution from electric and non-static magnetic gluons. The ratio &o/^o i^ ^ well-defined number and 
given in Eq. ^. These first two terms comprise the leading logarithmic approximation. The last term, 
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c, represents contributions that go beyond leading logarithmic accuracy, such as terms of order one in 
J d{q — uUeq, or of higher order, for instance ~ J d(q — fi). We have neglected these terms in our derivation 
of Eq. (^l]). In Eq. (||), they were written in terms of the undetermined constant 5q. We are therefore 
permitted to choose c = 21n6Q in the following. 

In order to understand which terms are responsible for the exponent and the prefactor in Eq. (nh let us 
solve ( [75| ) assuming that the gap function (j>q is constant as a function of momentum. This will give the 
wrong coefficient in the exponent, but suffices to understand the parametric dependence of (m on g. For 
k = n, (j)q = (f>k = 4>o = const., the {q — /i) integral in Eq. ( [75[) can be performed to yield {S ^ 0o) 



In" 



2S\ ,Jb^],J^S 



(76) 



The first term arises from the leading term ~ ln(/i^/|e^ ~ ^k\) in (uSl) , while the second comes from the 
leading logarithmic term ~ ln[6o/(&o.9^)]- The quadratic equation ( [Tq ) for \n(2S/(f>o) has the solution 



-I) 




(77) 



In weak coupling, the term ~ 2/^^ dominates the right-hand side, so that we can expand the square root. 
This term then gives rise to the exponent in 1/g for (j)o. The term ^ \n[bo^/{g^S)] in ( [77| ) gives rise to the 
prefactor of the exponential. In this way we obtain 



/.o=^2^A^expU^j x[l + 0(5)] . 



(78) 



Note that, to leading order in g, the dependence on the cutoff S cancels in the final result. Equation ( fTq ) 
is rather similar to Eq. (nh. The difference is that, due to our (erroneous) assumption of a constant gap 
function, the coefficient in the exponent is incorrect, V2, instead of 7r/2. 

B. Converting the integral equation into a differential equation 

The gap equation ( [7^ ) can be solved analytically by approximating the logarithm under the integral. As 
in M, we replace 



im 



,2, ,2 



6V' 



■-q ^k 

and the gap equation ( [72| ) becomes 
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bfi 



e/c / Jo 



^^(y) 9iq-k)+\n(^] 0{k-q) 



^ii^ tanh (±l\ 



\2TI 



d{q ~ n) 



k—}i 



tanh(^)lnr^ 



(79) 



(80) 



Upon differentiation with respect to k, we see that the gap function 0^ monotonously decreases from its 
maximum at the Fermi surface k — fi (we assume (fiq > for < q — fJ- < S) . Let us introduce the variables 



y = In 

a; = In 
X* = In 



2 6/i 



q- fi + eq 

2bfi 
k - fi + ek 



2bpi 



In 



2bfj, 



(81a) 
(81b) 
(81c) 



where (j) = 4>u = 0(e^, M^) is the value of the gap at the Fermi surface. Obviously, x* > x, y. As (5 <C TUg <^ 
IJ, -^ bfi [cf. (|73|)] , both X and y are much larger than 1. In terms of these new variables. 
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Cq — b^e 
ek = bfj.e~ 



l+fMe-i^'-y)\ 



1 



-(a:*— a:) 






(82a) 
(82b) 



where we defined a new function e(y). 

Tlie transformation to the new variable y is natural, because 



dy 



d{q - y) 



(83) 



is the measure for integration, without any further Jacobian. It is similar to Son's variable yson — ln(/i/(7o)- 
For Son's variable, however, dyson = —dqo/lo, so that in the gap equation, it is necessary to include a 
Jacobian for the transformation from qo to yson- This Jacobian does not affect his results for the parametric 
dependence of the gap function on g, but to leading logarithmic accuracy, it does affect the prefactor and 
the shape of the gap function, cf. discussion at the end of subsection III C . 
In the new variables x, y the approximate gap equation reads 



<t>{x) 



with e{y) given by 



:£ dy tanh (^^) 4>{y) + £^ ^ dy y tanh (^^ 



^{y) 



(84) 



(p2^. Here, we have neglected terms ~ In [l + {(j){y) / (j))^ e ^'■^ ^^] <C y and 
In [l + {(l){x) / (f))"^ e~^^^ ~^'] <^ X. This is a good approximation, as the value of the logarithms is bounded 



<j>ix) 



from above by In 2 ~ 0.693, while x, y > hi{bfi/6) 3> 1. Differentiating with respect to x yields 

d 
dx 

This is the generalization of Son's equation for the gap function to nonzero temperature. 



dy tanh ( ^ ) 0(y) , 



e{x) 



9 tanh — - 0(x) 



(85a) 
(85b) 



C. Zero temperature 



Th e ap proximation (79) has succeeded to convert the original integral equation into the differential equa- 
tion (85b). For nonzero T, however, because e{x) is a complicated function of x, this equation still requires 
a numerical solution. On the other hand, for T — the solution is simply a trigonometric function. Deter- 
mining its phase and amplitude from the values of (f>{x* ) and d(j)/dx at x* , we obtain 



(x) = 00 cos [g{x* - x)] 



(86) 



The value of the zero-temperature gap function at the Fermi surface, (j)o, can now be obtained by inserting 
the solution (pq) into the approximate gap equation for x = x* . This yields 



i>Q — 26 exp 



1 
-— arctan 
9 



1 



9 ln(&At/'5) 



(87) 



The dependence on the momentum cutoff d is spurious: in weak coupling, g <S^ 1, we may expand the arctan 
with the result 

00 = 2 6m exp f-^) X [1 + Oif)] . (88) 

This is identical to Eq. (m), except that here the undetermined constant fog is equal to one, since we computed 
only to leading logarithmic accuracy and dropped terms of order one in the gap equation. 

While our analysis is strictly valid only for small values of g, it is instructive to extrapolate to strong 
coupling. The behavior of </>o as function of g is shown in Fig. H^. One observes that 0o has a maximum at 
a coupling constant g ~ 4.2. The maximum value is quite large, 0o — 0.26//, due to the large prefactor bo, 
Eq. (0), which has important implications for phenomenology |15|. 
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0.20 




FIG. 1. <^o/m fi^s function of g. We have set the undetermined constant fog in Eq. (|2|) equal to 1. 
To leading order in g, x* = 7r/(2g), such that 

(l)(x) — 2bfj, exp I — -— ) sin(^a;) . (89) 

V 25-7 

The gap function peaks at the Fermi surface, a; = cc*, and varies over a region x* — x '^ I/5. This iniphes 
that the gap integral is dominated by momenta exponentially close to the Fermi surface, e^ ~ bfi e~^ <C mg. 
For example, when q — fi ^ rUg, and consequently y ~ 1, sin(gy) ~ gy, and the gap function is g times 
smaller than at the Fermi surface. Physically, this dependence of the gap function on momentum reflects 
the fact that single-gluon exchange is dominated by forward-angle scattering |13]. 



The effect of using the variable y, Eq. (|81a| ), instead of a variable z = ln(5/i/eg) analogous to Son's 
variable yson = ln(/^/go) [||i0 is the appearance of a prefactor 2 in Eq. (|8|), as well as a factor of 2 under 
the logarithm in x* = \n(2bn/(p). These factors of 2 were found empirically in a numerical analysis by Schafer 
and Wilczek [111. 



D. Nonzero temperature 

At nonzero temperature, the gap function depends on both x and T. In this subsection we denote this 
dependence by (/)(a;, T). The value of (j>{x, T) at the Fermi surface is denoted by 0(r). As before, 0o = 0(0). 

In weak coupling, we can compute 4>{x,T) by assuming 0(x,T) ~ (j){T) 4>{x,0)/(j)o- In other words, the x 
dependence of the gap function is taken to be the same as at zero temperature, so that the only effect of 
nonzero temperature is to change the overall magnitude of the gap function. We first present this calculation, 
and then discuss why, in weak coupling, it is reasonable to assume that the gap function as function of x does 
not change when the temperature is varied. Finally, we verify this assumption by numerical calculations. 

For reference, we review how the solution to a BCS-type gap equation changes with temperature. Taking 
the BCS coupling constant to be G, the gap equation is of the form 

</.(T) = G2^^fc^tanh(^)0(T) , e, = ./i^~;ir+Wy ■ (90) 

In BCS-type theories, the momentum dependence of the gap equation can be ignored, so that we just have 
a fixed-point equation which determines </>(T). Besides the trivial solution, ^(T) = 0, at sufficiently small 
temperatures there is also a non-trivial solution (f>{T) 7^ 0. It is determined by solving the equation 

l^G'£:^i:i-J^tanh{^) . (91) 

At zero temperature, the solution is 

l^G'lnf^] , </,o~2Aexp(-^) . (92) 
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The gap equation requires us to introduce a cutoff, A, but all that matters is that the cutoff is much larger 
than the value of the gap at zero temperature, A^ (pQ. 

With increasing T, the thermal factor tanh(eq/2T) reduces the integrand in (|9l|), so that (f>{T) has to 
decrease in order to achieve equality of right- and left-hand sides of Eq. (91). Above a critical temperature. 



Tc, this is no longer possible and we only have the trivial solution, (f>{T) — 0, T > Tc- 

Let us investigate in more detail how the balance between the left- and right-hand sides of Eq. ( |9l| ) is 
achieved at nonzero temperature. We divide the integration region in (vW into two pieces, for q — jjl smaller 
or larger than K(/)o, where k, is some pure number, which is assumed to be large. The first region corresponds 
to momenta near the Fermi surface, < q — fi < k^q. Although k 3> 1, because the gap is exponentially small 
in the coupling constant, Eq. (192), in weak coupling this region constitutes only a rather small contribution 



<a- 



< 



to the complete gap integral in (|9l|). 

The bulk of the integral comes from the region of momenta away from the Fermi surface, K(l 
A. The quintessential point of our argument, which we shall apply to the QCD case shortly, is that, for 
temperatures on the order of 0o, the thermal factor can be neglected in this region, up to corrections 
^ exp{—K(j)Q/T). The gap equation (pi|) can therefore be written in the form: 



1~G^ 



K<i>0 



d{q - fi) 



"■"•(I?) 



d{q - fj.) 



K0O 



(93) 



In the small region very near the Fermi surface, for momenta < q — fi < k(J)q, the thermal factor must be 
retained. 

In weak coupling, G <C 1, the region away from the Fermi surface almost saturates the 1 on the left-hand 
side of the gap equation (pq). To see this, we compute the contribution of this region to the right-hand side. 



G 



2/ ^(i^^G^i^M )^i_G^in(2.) 



(94) 



where we used the solution ( |92| ) at T = 0. Hence, in order to fulfill the gap equation (|93D, the first integral 
in (93) has to balance the term G^ ln(2K), 



G'^ ln(2K) ~ G^ 



d{q - fi) 



tanh 



'0 tg 

This equation can be written in a more concise form. Note that 



V2r 



^ 



d{q - ^J.) 



Iii{2k) 



(95) 



(96) 



where e" — ^/(<7 — /j)^ + (pn is the quasiparticle excitation energy at zero temperature. We therefore replace 
the term ln(2K) in (Bsh by this integral. 



0~G' 



K0O 



d{q - fi) 



tanh 



(1^\ 



(97) 



The same result could have been obtained directly, by subtracting Eq. (91) at zero temperature from ( |9l| ) at 
T 7^ 0, and using the fact that, for momenta g — /z > Kipo, eq — e°, so that the contributions to the integrals 



away from the Fermi surface cancel in the subtraction process. The integral in (97) is finite in the infrared, 
even for (l>{T) — 0, since tanh(x)/x ^ 1 as x ^ 0. 

Dividing all quantities with dimension of energy by (po, one realizes that Eq. ( |97| ) determines the ratio 
<j>{T)/(j)o as a function of T/4>o. Note that this ratio is independent of the cutoff A, as well as the coupling 
constant G, since G^ is just an overall constant in (p^) which we can divide out. 

Equatio n (|97| ) has the following interpretation. As mentioned above, the 1 on the left-hand side of the gap 
equation (|93|)is almost completely saturated by momenta away from the Fermi surface, kc^q < q — h < A, 
where thermal corrections are negligible. Thermal effects become important in a small region near the Fermi 
surface, < g — /i < K(j)o. Only in this region, the change of the gap function with temperature has to 
compensate for the presence of the thermal factor. The ratio 0(T)/0o is therefore completely determined by 
investigating how the gap equation changes with T in a small region around the Fermi surface ! 
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In Eq. (p7|), we are allowed to send k — > oo, as done in Eq. (10) of ||l^, because tanhx -^ 1 ior x ^>- oo, 
and Eg — > eH^for g — /i — > cx). The ratio 0(r)/0o is therefore not only independent of the cutoff A, and the 
coupling constant G, but also of k. However, sending k — > oo somewhat obscures the fact that only a small 
region around the Fermi surface controls how (f>{T) changes with temperature. 

Equation ( |97| ) cannot be evaluated analytically for arbitrary T, although it is easy to solve numerically; 
as is well known, the transition is of second order, with 4>{T) — > as T — > T^. At the critical temperature 



Tc, however, one can evaluate (|9^) analytically: 
ci(g-^) tanhf^]-- 



Here 7 ~ 0.577 is the Euler-Mascheroni constant. The solution is 

11^^^ 0.567 . (99) 

00 2 

which is the usual result in BCS-theory ||^,g] . 

We now go through a similar computation for QCD, neglecting the change in the x dependence of the gap 
function with temperature, cj){x,T) ~ (j){T) (j){x , 0) / (f>o ■ Consider the gap equation ( p4[ ) for x — x* , i.e., for 
the gap function at the Fermi surface, (/>(T). Assuming T < T^, and taking the non-trivial solution (f>{T) 7^ 0, 
we can divide both sides by (j){T). We obtain 

-2 r .,.„._.^(y)^ 0(y,o) 



1 ~ 5M dyy tanh ^ ^^^^ . (100) 

This is the equation analogous to (pi|). We now assume that the critical temperature is on the order of the 
value of the condensate at zero temperature, 0o- In analogy to the treatment for BCS-like theories, we divide 
the region of integration into one in which < g — /i < k(J)o, and one in which kc/jq < q — ^ < 5. As before, 
for momenta away from the Fermi surface the thermal factor is neglected. Using (j)[y , 0) / (/)q = sin(gy), we 
obtain 



f 



/ dyy tanh ( — - I sin{gy) + dyy sm{gy) 

Jx^ \ ^J- J J\n{bu/S] 



(101) 



which is analogous to (|93[). Following the definition of the variable x in ( 81b ), here we introduced x^ = 
ln[6/i/(K(/)o)]. Again, in weak coupling the 1 on the left-hand side is almost completely saturated by the 
integration region away from the Fermi surface. To see this, we compute the respective integral in a power 
series expansion in g: 

f dyy sm{gy) ~ 1 - - g H2k) + 0{f) . (102) 

J\n(bii/S) ^ 

In order to fulfill ( tlOlf), th e integral over the region very near the Fermi surface has to compensate the terms 
of order 0{g) in Eq. (jl02|). 



f j^ dyy tanh (^^ j ^xn{gy) ~ | g \n{2K) + 0{f) 



(103) 



To see how this compensation works, we expand the various terms on the left-hand side in powers of g. As 
X* = TT / {2g), and x,^ = x* — ln(2K), y is of order a;* in the whole integration region, up to corrections of order 
g. Note that this is equivalent to approximating the factor ln(6/x/eg) by ln(26/x/(/)o). For the momentum 
dependence of the gap function, this has the consequence that sin(gy) = 1 + 0{g), i.e., to leading order in 
weak coupling, the gap function can be taken to be constant, (j){y, 0) ~ (/>o. In this way one obtains 

f -9 ln(2.) . I , f' to^ tanh (g.) + Oif) , (104) 

where we reverted the integration variable y into q — /i. With (|96|), we can write this as 

18 



0{f). (105) 



Apart from the prefactor (which can be divided out), to leading order in weak couphng, Eq. ( |105| ) is the 
same as in BCS-type theories, Eq. (£^) ! Thus, also the ratio (J){T)/(J)q is identical to that in BCS. Given the 
very different nature of the gap equation, this is a remarkable result. 

How can we claim that we can reliably compute (j){T)/cj)Q to leading order in weak coupling, although we 
cannot compute the overall magnitude of the condensate at zero temperature, the constant 6q in (||) ? To 
understand this, consider the coun ting of powers of g in the gap equation; for this, it suffices to use the 



approximations of subsection [II A 



From Eq. (uw we have seen that the exponential in 1/g arises from terms ^ 3^ In {2S/(j)o). Since 
3 ~ exp(— 1/5), these terms are of order 1, and therefore of the same order as the left-hand side of 
(yW- Analogously, in the nonzero-temperature gap equation (101), terms of order 1 arise from the region of 



lich 



momenta away from the Fermi surface, and thus balance the 1 on the left-hand side. As temperature effects 
are negligible in this region, we conclude that the exponential behavior '^ exp(— 1/g) of the gap cannot 
change with temperature. 

The prefactor of the exponential is determined by terms ~ g^ \n(26/(/)o) ^ g in the gap equation 
These include all terms of t his o rder, i.e., the leading logarithmic terms as well as terms of order one, wr 
give rise to 6q. From Eq. (|103| ) we realize that temperature effects also enter at order g. These effects 
therefore change the prefactor of the exponential. In the above discussion, we have determined the change of 
the gap equation with temperature to leading order in g, or in other words, we have determined the prefactor 
at nonzero temperature. To leading order in g, the result is that the prefactor changes precisely in the same 
way as in BCS-like theories. 

To understand how the gap function changes with both x and T, consider first the region just below the 
critical temperature. Even though the overall magnitude of the gap, (f>{T), is small, we can still consider 
4>{x, T)/(j){T); even as T — *■ Tc, this ratio remains of order one. Then near the Fermi surface, the gap function 
must change due to thermal effects: after all, the variable x*{T) = ln[26/i/0(T)] diverges as T ^ T^ when 

To understand the change in the gap function, consider energies small relative to the temperature, e(x) < 



T. In this limit, the thermal factor tanh[e/(2T)] ~ e/(2T). Using the definition of the variable x, Eq. (85b) 
is approximately 

ej^ cbie,T)^-f[^) cj^{e,T). (106) 

About small e, the solution is given as a power series in e/(2T): 

He,T) ^ ^l-fJ- + 0{g*)\ct>{T) , T ^ T, , e«T. (107) 

Thus when the energy is much less than the temperature, although cc ^- 00, the gap function is not ~ sm{gx), 
as for the zero-temperature gap function, but approaches a constant. That is, due to the thermal factor in 
the equation for the gap function, the momentum dependence is cut off, and the gap function "flattens" . 
This is why in our previous calculation, we could neglect the momentum dependence of the gap function 
when the energy is less than k^q. 

Our analytical results are confirmed by numerical solutions of Eq. (B5q) at nonzero temperature. In Fig. 
y we show the gap function at the Fermi surface as a function of temperature. We considered two values of 
the coupling constant, g — 0.1 and 4.2. The first is safely in the weak-coupling regime, while the latter is the 
value where 0o has a maximum as a function of g, cf. Fig. |l|. In weak coupling, (j>{T)/(f)Q is indistinguishable 
from BCS. Surprisingly, even for large g — 4.2, while Tc is slightly larger than in BCS theory, Tc/4>o — 0.585, 
(j){T)/ (j)o is rather similar to the behavior in BCS-like theories. This is not unlike the situation in strongly 
coupled BCS, i.e. Eliashberg theory pj, where Tc/4>o also changes slightly, albeit in the opposite direction. 

In Fig. 11 (a) we show (l>{x, T) as a function oi x aX, g = 4.2 for two temperatures, T = and T = 0.581 (J)q 
which is close to Tc- As the temperature increases, the overall scale of 4){x,T) decreases, because the 
condensate is evaporating, cf. Fig. ^. Simultaneously, x*{T) = ln[26/i/(/)(r)] increases. In Fig. ^ (b) we plot 
4>{x ,T) / (j){T) as a function of x/x*{T) at the same two temperatures. As can be seen, once we divide out 
the overall scale 4>{T), the ratio (p{x,T)/(j){T) is relatively insensitive to temperature. At T = 0.58100: we 
observe that the gap function does "flatten" as it approaches the Fermi surface, x -^ x*{T), as we argued 



previously on the basis of Eq. (107) 
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FIG. 2. The gap at the Fermi surface 0(T), normalized to its zero-temperature value (/)o, as a function of T/^o- 
Solid: p = 0.1, dashed: g = 4.2. 
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FIG. 3. (a) The gap function (j){x,T), normalized to the zero-temperature gap function (j)o, as a function of the 
variable x for T = (solid) and T = 0.581 (po (dashed), g — 4.2. (b) The same, but now the gap function is normalized 
to its value <j){T) at the Fermi surface at the same temperature, and plotted as function of x normalized to x*{T). 
For the sake of simplicity, we have chosen 5 = 6/i. 



IV. CONDENSATES WITH TOTAL SPIN J = 1 
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A. Classification 

The gap matrix $+ for a condensate with total spin J — has the form (|l^) . Generahzing this equation to 
a condensate with total spin J = 1 is straightforward. For a, J — 1 condensate, the individual gap functions 
0^^ in Eq. (H) are no longer scalar functions in coordinate space, but 3- vectors, 0^^^. These vectors have a 
component along the direction of k, 

ipUK)^4>lsiK)-k, (108) 

and two components transverse to this direction, 

xLW^0LW-(l-kk), (109) 

such that 

0LW^^LWk + xLW- (110) 

The gap matrix <&+ is, however, still a scalar in coordinate space. One therefore has to contract the vector 
indices of Xhs with the other independent 3- vector at our disposal, 7, 

xlsiK) ■ 7 ^ cPUK) . (1 - kk) . 7 ^ xL(^) • 7±(k) , (HI) 

where we have defined 7_L(k) = (1 — kk) • 7 and used the projector property (1 — kk)^ = 1 — kk. 

Suppressing color and flavor indices for the moment, the most general ansatz for a gap matrix describing 
J = 1 condensates is therefore 

J = 1 : <i>+{K) = E E E i^'^siK) + xLiK) ■ 7±(k)] V^O^ . (112) 

h=r,i s=± e=± 

The scalar product xtsi-^) ' 7-L(k) has simple commutation properties with the quasiprojectors T'^g(k): 

Vh7±{k)='r^{k)V-h. , (113a) 

P,(k)7i(k) =7i(k)P_,(k) , (113b) 

A^(k) 7i(k) = 7i(k) [A^(k) - e ak7o] ■ (113c) 

In the massless limit, a^ = 0, and 7j_ commutes with the energy projectors. 



Instead of contracting the transverse components (109) of (f)'^^{K) with 7, one could have alternatively 



also used (^^^^{K) itself. The commutation properties for cpf^^ ■ 7, however, are more complicated than for 
Xhs ' 7-1- (k), with additional terms proportional to k. The advantage of using 7±(k) is that these terms are 
cancelled by the projector 1 — kk. 

Finally, an important question that arises is why we choose 7 to contract the transverse components of 
(/>^g in Eq. ( |112| ). Why not 707, or 75707 ? The answer lies again in the com mutation properties ( |113| ). 
In the commutation property for 70 7±(k) or 75 70 7_L(k) corresponding to ( 113c| ), the sign of the energy on 



the right-hand side would be flipped. Physically, this means that such ter ms r epresent pairi ng of p articles 



with antiparticles, and not of particles with particles, cf. discussion of Eq. ( |119| ) in subsection IVC. This is 
clearly not what happens in a superconductor, where particles form pairs near the Fermi surface. 

B. Parity 

The J = 1 condensates can be decomposed according to their transformation properties under parity. Let 
us first write the interaction term ipc $'*' "0 in the effective action |9[| in the form 

MK)<^+iK)^{K) = E E E MK)cI>UK) ■ (k + 7±(k)) n.(k)V'(i^) . (114) 

h=r,e s=± e=± 

Let us define new condensate fields 
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0+ ^ i (0++ ± 0+_) , 0± ^ ^ (c/.;_ ± 07+) 



7r+ ^ 1 {<!>+_ ± 4>t 



In terms of these fields, the interaction term (114) assumes the form 



e=± 



J2^ciK) \ </.;(i^)-(k + 7i(k))+01(i^).(k + 7^(k),75 
n%{K) ■ (k + 7+(k)) + n'L{K) ■ (k + 7i(k)) 7 



7r± ^ - (0-+ ± 07_) . (115) 



1 + 


6 70 7 


k 


1- 


2 

6 70 7 


k 



A^(k) 



A'=(k) ^{K) . (116) 



In the Hmit ?ti ^ 0, A'^(k) ^- (1 + 6707 • k)/2, such that the terms proportional to tt^ vanish. 

Under a parity transformation, 70 ^ 7o, 7 — > —7, 75 — > —75, and k — > — k. Thus, the terms (lie 70 7-k)/2 
and the energy projectors do not change under a parity transformation. The parity of V'c is opposite to that 
of ■0 g. A parity transformation leaves the effective action invariant, since the QCD Lagrangian (from which 
the effective action is der ived ) does not violate parity explicitly. Thus, if we perform a parity transformation 



of the interaction term (116), the gap functions 0j_, tt^. must have the same parity as the accompanying 



Dirac matrices. Consequently, we conclude that cj)'^, tt^ are parity-even, while (p'L, ttI are parity-odd. 

Two limiting cases are of special interest. If all right-handed, positive- helicity gaps are equal to the left- 
handed, negative-helicity gaps, (pf.^ = 4>£_, and all right-handed, negative-helicity gaps are equal to the 
left-handed, positive-helicity gaps, 0^_ — 4'1+^ then 01 = ttI = 0, i.e., the odd-parity gaps vanish, and 
condensation occurs exclusively in the J^ = 1+ channel. [It is proven in Appendix |A| that in this case Eq. 
(116) agrees with the ansatz of Bailin and Love for the J^ = 1+ gaps B.] On the other hand, if they are 
equal in magnitude, but different in sign, condensation occurs exclusively in the J^ = 1~ channel. 

This is different for the J = gaps. In that case, one can write down an equation analogous to (116). 



The difference is that all vectors 0j_, tt^. are replaced by scalar functions 0j_, ttj., and the terms k + 7±(k) 
are absent, too. For the parity transformation properties, this has the consequence that if 



±0|, , condensation occurs in the J^ — 0^ channel. 



'r+ 



±(, 



C. The ultrarelativistic limit 



In the following, we exclusively consider massless fermions. In this case, the helicity projectors (and 
corresponding indices) become superfluous, cf. Eq. (|25|), 



TO = 



^hs 



vl 



Xhs 



Xh 



(117) 



and the quasiprojectors become true projectors, cf. discussion above [and Eq. (B29) ^\. The general ansatz 
(112) simplifies to 



J = 1, TO = : ^+{K) = E E ['^^ W + XI{K) ■ 7^(k)] P;j(k) 



(118) 



h=r,t e=± 



In the effective action, the interaction term V-'c ^^ V' decomposes into four parts on account of the projectors 
VI [cf. Eq. (B34) of §], 

MK)<^+{K)^{K)=Y, Y,[rcf,{K)^l{K)rH{K)+rc-h{K)xl{K)-^^{)^)MK)\ . (119) 



Here we defined V'/j = ^^ V': and used the commutation properties ( |113| ) of 7^. From (119) it is obvious 
that the longitudinal gaps correspond to condensation of quarks with the same chirality, {rr) or [U). In 
this respect the longitudinal J = 1 gaps are similar to the J = gaps [g|. On the other hand, the transverse 
J = 1 gaps c orre spond to condensation of quarks with different chiralities, {rt} or {Ir). 

Equa tion ( 119 ) shows why we do not choose 707 or 75707 to contract the transverse components of 0^^^ 



112|), cf. discussion at the end of subsection IV A . The extra factors 70 or 75 70 flip the sign of the energy 
113q) , and the transverse gaps would describe pairing of particles with antiparticles of the same chirality. 
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D. Color and flavor representations 

With the help of the symmetry property [cf. Eq. (B4) of M] 



C<^+{K)C-^ = [^+i-K)]^ , (120) 



[Vl{-K)f = -'PI(K) , [xl(-K)f = -X-hiK) ■ (121) 



we derive 



The symmetry properties (t21j) allow us to classify the possible color and flavor representations of the J = 1 
condensates for massless fermions. We assume there are Nf flavors of massless fermions with a global flavor 
symmetry SU{Nf)r x SU{Nf)i and, of course, a local color symmetry SU{3)c- From group theory, 

2x2 = 1^ + 3, , 3x3 = 3a + 6, , (122) 

where the subscripts a and s denote antisymmetric or symmetric representations, respectively. For single- 
gluon exchange, the color-antitriplet channel 3^ is attractive, and the color-sextet channel 6^ is repulsive. 
We first discuss the longitudinal condensates (pf^. From Eq. ( |ll9| ), the longitudinal J = 1 condensates only 



couple quarks of the same chirality, and so do not break the SU{Nf)r x SU{Nf)i fiavor symmetry. By (121) 
the longitudinal condensates must correspond to a color-flavor representation which is overall antisymmetric. 

For Nf = 1, the flavor representation is trivial, and condensation must occur in the 3^ channel. This 
is most likely the favored channel for condensation of quarks of the same fiavor, since the J = gaps are 
overall symmetric MM , and consequently must be in the repulsive 6^ channel. 

For Nf = 2, the allowed color-flavor representations are either (3^, 3{), which is favored, or (6^, l{). This 
is in contrast to the J = gaps, which come in (3J^, 1^) or (6^, 3{). Here the fiavor representation refers to 
either SU{2)r or SU{2)e. _ 

For Nf = 3, there are (3^^, 6{) or (6^, 3^). On the other hand, the J = gaps are (3Jj, 3^) or (6^, 6{). 

If only longitudinal J = 1 gaps are present, the parity properties are analogous to the J = gaps. In 
particular, since condensates of different chiralities do not mix, the magnitude of the longitudinal J = 1 gaps 
will be equal, while their relative phase represents the spontaneous breaking of parity [7|Jl4[| 



The symmetry relation (121) relates transverse J = 1 condensates of different chirality. Therefore, we 
cannot draw general conclusions about the possible color-flavor representations of the transverse J = 1 gaps. 
A priori, both the symmetric as well as the antisymmetric color and fiavor representations are allowed. 

Two special cases are of interest. H Xr i = Xf r' ^*^ ^^^^ condensation is in the J^ = 1+ channel, the 
color-flavor representations of the transverse J = 1 gaps are identical to those of the longitudinal J = 1 
gaps. On the other hand, ii Xr e= ~X'e r' s° ^^^^ condensation is in the J^ = 1~ channel, they are equal to 
those of the J = gaps. 

When Nf > 2, the appearance of transverse J = 1 gaps must necessarily break the SU{Nf)r x SU{Nf)i 
symmetry to a vector-like SU{Nf) symmetry. This striking feature arises because the transverse gaps are 
proportional to 7^, and not, say, 70 Tx- 

E. The gap equation 

In general, condensation can occur in channels with arbitrary tot al sp in J . Therefore, the gap matrix 
^+(K) wiU not simply be of the J = form (ph, or the J = 1 form (pT2|), but it wifi be a sum of (fgh and 



(112), as well as contain contributions from higher spin J > 2. We do not attempt to solve this problem in 
full generality. Instead, we consider the simpler case where the gap matrix contains just J = and J = 1 
contributions. We also take the fermions to be massless. The gap matrix then assumes the form 

J = and 1, m = : ^+{K) = J2 J2 [<^'^(^) + ^hiK) + xl{K) ■ 7^(k)] V^O^) . (123) 

h=r,t e=± 

where (f>\ denote the J = gaps, and Lp\ the longitudinal and xt, '^^^ transverse components of the J = 1 
gaps. 

The quasiparticle propagator can be computed from Eqs. (R), (H), and (S), and the commutation property 



(113c). Ignoring the color and flavor structure for the moment, one obtains 

23 



G+(Q)-<( E El'^o-l'Z-eM)' (124) 

h=r,l e=± 

+ ( [</.^(Q)]t + [^^,(Q)]t - [xl,(g)]t • 7±(q)) ( ^1{Q) + ^liQ) + xliQ) ■ 7±(q))] ^(q)}"' [G^iQ)] ' ■ 

To invert the term in curly brackets, we first have to specify the color-flavor structure of the various con- 
densates. In the following, we just consider the most simple case Nf = 1. Let us assume that condensation 
occurs exclusively in the attractive color-antitriplet channel 3^, i.e., <i>^- = ^ijk^X ^^ ^^^ antisymmetric 
Nc X Nc matrix in color space, i, j = 1, • ■ ■ , Nc- In this case, the J = gaps vanish, 
of Fermi statistics 



= 0, on account 
Since the individual gap functions (/7| ^ ■ , Xh u ^^^ ^^^o in the 3^ representation, we 
conclude from (121) that x^ = yf. I n contrast, ^'^ and (^| remain a priori unrelated. 

The term in curly brackets in (124) has an off-diagonal contribution in color space. Since its contribution to 
G^ is quadratic in 0^, we shall neglect it in the following, and consider only the diagonal part, G^- ~ 5ij G+. 
Even so, the inversion of the term in curly brackets is still cumbersome, due to the presence of terms ^ 'y ±^. 
To simplify the treatment, we assume that the J = 1 gaps are real-valued. With Eq. (121) this leads to 



[^1{K)]^ ^ Wli-K)]" ^ ~^liK) , [x'_,iK)]^ = [x'-d-K)]"^ ^ -XI{K) 



(125) 



This has the consequence that all cross terms ~ 7_l between the gaps 93^ and Xh vanish in (124). The 
quasiparticle propagator can now be explicitly computed: 



G+(Q)-E E 



vfS^ 



h=r,i e=± Qq 



m)]' 



[GoiQ)]'" 



(126) 



The quasiparticle energies are 



<(0I) ^ 



Na 



■'hi 



1/2 



(g-eM)2+^0|,.c/,; 
/=i 

Inserting (126) into the gap equation (g), the same steps that led to ( |2q ) now lead to 



^^|E7'^A,.(i.-Q)E E ^"^^l"^".^'^^l'7"^^^ ^-i(q)7- 

Q h=rj e=± 
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(128) 



The gap equations for different fundamental colors i decouple, therefore we omit the color index in the 
following. Taking projections we arrive at the following equation for the longitudinal gap functions: 



Q I 



vliQ) 



<l'o 
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q'o - [^^'(.<pr)] 



Tr[7';;(k)7^7'ir,(q)7l 
-Tr[7';^(k)7^7'!:,,(q)7l 



(129) 



This equation is rather similar to Eq. (M) . The only difference is th e appearance of the transverse gaps in 
the quasiparticle and quasi-antiparticle energies e^((/)|), cf. Eq. (127). 
For the transverse gaps, we use the fact that xt = Xe to arrive at 



xliK) 



2 oT 

- 9 — 
3^ V 



EA^.(i^-g)| — 

Q [% 



xUQ) 



+ 



[e?('/'^J]' 



.Tr[7A-^(q)7"^A^(k)7^(k)7^] 
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(130) 



Taking Coulomb gauge for the gluon propagator, Eq. (|30|), and computing the traces similar to Eqs. (|3l|), 
we obtain 
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^liK)^lg'^Y. 
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^?(0I)]^ 
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qi-bn^H^)]' 



^ 1 + k ■ q ^^ k + q 



Ai{K - Q) - At{K - Q) 1- 
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(k-q)^ 
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(131) 



Neglecting the antiparticle contribution, and taking fc = g = /i, as well as k ~ q (the collinear limit; in Eq. 
( [70| ) this corresponds to neglecting terms of order p^), we obtain for the transverse quasiparticle gap 



xtiK)^l9%j: 



xtiQ) 



<ll-[4i^t 



[A,{K - Q) ~ At{K ~ Q)] 



(132) 



To leading logarithmic order, the gap equation for the tra nsve rse gap is identical to that for the longitudinal 
gap. To see this directly, one computes the traces in Eq. (12£). One obtains an equation very similar to Eq. 
(p3). Now one computes the angular factors for k ~ q, k — q = ^■. 



1 



3-k-q 1 + k-q (fc-q)2 



(k-q)2 



(133) 



which proves our assertion. 

Consequently, to this order all J = 1 gaps are of equal magnitude. Further, the sum of the squares of the 
J — 1 gaps fulfills the same gap equation as the square of the J — gap. This conclusion agrees with Son's 
renormalization-group analysis ^ , who argued that the parametric dependence on g of any spin J gap is 
the same. We find that even the prefactor is the same to leading logarithmic accuracy. Our results differ 
from those of Hsu and Schwetz |1^ , who argued that J = gaps are favored over those for higher spin. 

Beyond leading logarithmic order we suggest that, if J = gaps are allowed by color- flavor symmetry, 
they are probably favored over the J = 1 gaps. If only a J = 1 gap is allowed, as for Nf = 1, we believe 
that either the longitudinal gaps, or the transverse gaps with a definite color-flavor representation, will be 
favored. Which one is favored will be determined by the constants 69 and b'^ in Eqs. (g) and (||). 



V. CONCLUSIONS 



We conclude by addressing effects which can contribute to the constants 6q and b[ in the condensate, Eqs. 
W) and (||). We suggest that our lengthy calculations which led to the result (0) for bo may be done much 
more efficiently by constructing an effective theory for quarks near the Fermi surface, as initiated by Hong 
pO| . This is presumably the easiest way to calculate 69 and b'l as well. Nevertheless, we can estimate what 
kind of effects could contribute to these as of yet undetermined constants. These include: 

(i) Gauge dependence of the condensate: 
If calculated properly, any physical quantity must be independent of the choice of gauge. For color super- 
conductivity, what is physical is the gap on the quasiparticle mass shell. At nonzero temperature and zero 
quark density, general arguments due to Kobes, Kunstatter, and Rebhan |2^] indicate that the mass shells 
for quarks and gluons are gauge invariant. Their proof does not extend obviously to nonzero density, but we 
shall assume that to be the case. 

At higher orders, the quasiparticle self energy, S+, contains not only the interaction with the condensate, 
Eq. (p), but also wave- function renormalization. Any apparent gauge dependence in the former must be 
cancelled by the gauge dependence of wave-function renormalization, leaving the quasiparticle mass shell 
gauge invariant. 

It is easy to see that such subtleties do not enter at the order to which we have computed. The gluon 
propagator in Coulomb gauge, with gauge fixing parameter ^c is given in (pQ). Alternatively, one could 
have taken covariant gauge; with gauge fixing parameter ^, the covariant gauge propagator was given in Eq. 
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(37) of [2^. We have seen, however, that to the accuracy to which we compute, only the gluon propagator 
in the static Umit, po — > 0, matters. In the static hmit, the gauge dependent terms are identical for either 
Coulomb or covariant gauges, and appear only as spatially longitudinal terms, ~ ^P^p' /p^- Consider, 
for example, how the gauge dependent part of the gluon propagator affects the quasiparticle contribution, 

~ (j>h/ \qx)~ \_^qi.(t>h)\ f' to the gap equation for (j)'j^. From the gap equation (|2|), and using (^, this 
becomes 

ic Tr [P+(k) 7 • PPI, (q) 7 ' p] ^ -& ^-^^ ^^^ - , k,q^^,. (134) 

z p 

That is, there are gauge dependent terms, but they only contribute to the antiparticle gaps, (jif^. (Further, 
the antiparticle gaps must be computed on their proper mass shell. At the Fermi surface, e^ ~ 2 /i is not 
small, in contrast to e^ ~ 0o ^ A*-) Consequently, gauge dependent terms in the particle gaps do not appear 
to even affect the prefactor in the condensate, the constant h'^. These conclusions about gauge invariance 
agree with the results of Schafer and Wilczek Q. In contrast. Hong et al. IQ argue that Landau gauge is 
preferred, as in other approximate treatments of Schwinger-Dyson equations. We insist that in the present 
example, direct calculation demonstrates gauge invariance without further ado. 

(ii) Wave-function and vertex renormalization: 
As noted by Son |q], one can have infrared singular factors for wave- function renormalization. For non- 
relativistic fermions, this was noted long ago by Holstein, Norton, and Pincus |2^. (This wave-function 
renormalization is not the HDL correction discussed by Schafer and Wilczek |13|; such corrections involve 
two hard lines, and are down by g^.) The dominant corrections involve a very soft transverse gluon on a 
quark line; this produces gauge-invariant terms of the form Z — 1 ^ g^\a[^/eq) ~ g when tq ~ (J3a- This 



correction was computed by Brown, Liu, and Ren |17|, who find that it is a large effect. Besides such wave- 



function renormalization, one might expect that the Slavnov-Taylor identities would also generate similar 
corrections for the gluon, and for the quark-quark-gluon vertex. This was not found, however, by the authors 
of0. 

(in) Effects of the condensate: 
We have computed the gap using an HDL-resummed gluon propagator. This is possible because the momenta 
which generate the gap are much larger than the scale of the gap. To see this, note that in the transverse gluon 
propagator, (p8[). Landau damping contributes to the gluon propagator when the momentum p^ ^ {m^Lo)'^\ 

since the frequency uj ^ (f), the dominant momenta are p ~ nig (f>^'^ , which for small (/) is much larger than 

Effects of the condensate on the gluon propagator can be estimated by power counting at large momentum, 
p 3> (po- One would naturally expect that they are ~ 5^'/'^, but due to an infrared divergence, they are larger, 
~ m'i(f)/p |20[| . These terms are important when m'i4)/p ~ p^, orp ~ mg (p^'^. This is exactly the same scale 
at which Landau damping operates. These effects will not alter the coefficient of the logarithmic divergence 
(and hence the exponent), but they will produce terms of order one in the gap equation, which contribute 
to 6q and b[. 

(iv) Damping of the condensate: 
In the above, we neglected the fact that the gap function has an imaginary part. To understand this 
imaginary part, consider first the self energy for a quark in a Fermi sea. As computed by Le Bellac and 
Manuel, and by Vanderheyden and OUitrault ^6| , away from the edge of the Fermi sea, the quark can decay 
into another quark and a very soft gluon. This is only possible with a HDL-resummed gluon, whose spectral 
representation has support from Landau damping in the space-like region. The damping rate of the quark 
behaves as r^ g^\p — ^\, vanishing at the Fermi surface. 

From a similar physical process, the gap function acquires a nonzero imaginary part when its momentum 
is away from the Fermi surface. A quark can scatter into a quark with a different momentum through a 
very soft gluon. We can estimate the resulting imaginary part of the gap function as follows ||lq| . If we had 
not restricted our analysis to the principal value part of the energy denominators arising in ( |42| ) , instead of 
In \e^ — e^l in Eq. ( |7^ ) we would have obtained ln(e^ — e^.). This logarithm has a cut for e^ < e^, giving rise 
to an imaginary part for cpk. 



-2 / rfeg , _ -2 , / efe 



lmcf>k--fl''^cf>g~f\n[^]cPo. (135) 
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Taking Ck ^ hfie ^, cf. Eq. (81b), momenta exponentially close to the Fermi surface occur when x 



t: / {2g) (in weak coupling). In this region, the imaginary part of the gap function, Im^fc ~ g {x* — x) ^o, is 
down by g relative to the real part, Rei/)^ ^ sin(^a;) (j)o. At the Fermi surface itself, x — x* ^ the imaginary 
part vanishes. Away from the Fermi surface, e^j. ~ /i, so a; ~ 1, and 4>k is strongly damped, with the real and 
imaginary parts of comparable magnitude, Ret/)^ ^ InK/)^ ^ g(j)Q. 

A gap function with a nonzero imaginary part is actually well known from strongly coupled supercon- 
ductors, as studied in Eliashberg theory [p|. Damping occurs for a similar reason as here, due to a nonzero 
imaginary part for the plasmon. 

There is no problem in principle with including the damping of the gap function. A spectral representation 
for the gap function is introduced, analogous to that of the quark and gluon propagators. The resulting gap 
integrals are more involved (especially at nonzero temperature), but can be treated in the manner which we 
employed above. 

(v) Magnetic mass: 
As argued in Section I, at zero temperature the scale for the magnetic mass is ~ /x exp(— l/g^). It is therefore 
negligible compared to the scale of the condensate. At nonzero temperature, the scale is no larger than g^T. 
For T r^ (f>Q, this is down by g^, and will only affect the prefactor of the gap to higher order in g. 

We conclude by stressing that the determination of the prefactor is not merely an interesting problem 
in its own right, but because it truly determines the physics of color superconductivity (at least in weak 
coupling). To the order at which we compute, there is absolutely no preference for the condensate to favor 
spin-zero over spin-one (or spin- two, etc. !). Surely the spin-one condensate is less favored than spin-zero; 
the ratio of the two condensates is, in weak coupling, a pure number which can be uniquely computed, once 
one knows how to compute the prefactor in the condensate. 

Indeed, perhaps one should entertain a more speculative hypothesis. Even if a J = gap is favored, maybe 
there is always some small admixture of higher-spin gaps, and rotational invariance is inevitably broken in 
the true ground state of color superconductivity. 
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APPENDIX A: THE J^ = 1+ GAPS 



According to the results of subsection [VB, if condensation occurs exclusively in the J^ = l"*" channel. 



the ansatz for the gap matrix reads (we suppress the dependence of the gap functions on K in the following) 



j^ = r 



$^ 



E 

e=± 



</>+ • (^k + 7_L(k)J + 7r+ • (^k + 7_L(k) 



A^(k) . (Al) 



Now use 



(k + 7.(k))i^^|^A^(k)=ki^ 



^4 



(l±/3kTak)-(l±/?k±ak)7-k 



(1 ± /3k) (1 ± e 7o 7 • k) =F ak 7 • k + e Qfk 7o 



(A2) 



to obtain 



$^ 



1 + /3k - Ckk ^^ + 



m + 0^ 



1 - /?k + ak 



(77+ + TT^ 
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) 


•k7o 


4 (<^+ + <^+) 4 (^+ + ^+) 


•k7k 


l + /3k + ak.,+ ,-x 1 " /3k - "k . + -X 
4 m ^+) 4 (^+ ^+) 


•k7o7-k 


^ (0+ + 0+)+ ^ (7r++7r+) 


•7 


4 W </'++^| 7r+)-77o 


^(-0+-0;+7r++7r;).77-k 




• 7 70 7 • k . 



Bailin and Love's ansatz for a J^ = 1+ gap reads [cf. Eq. (4.24) of 
dependence of the gap functions] 



(A3) 



we again suppress the momentum 



Ai • 7 + A2 • k 7 • k + i A3 X k • 7 75 + A4 • k 

A5 -770 + ^6 ■k7 •k7o + A7 • k7o + i As X k-77o75 



(A4) 



[We added an i in the last term as compared to [Q. This simphfies the notation in the foUowing, but is not 
essential, as the gap functions arc in general complex-valued.] With the definition of 75 — 1707^7^7'^ one 
computes 



i A X k • 7 75 



7o(A-k + A-77-k 



-70 A-7^(k)7-k 



(A5) 



and rewrites ( |A4| ) as 

A 



(A4 



(A7 



k7o + A2-k7-k-A6-k7o7-k 



Ai -7 + A5 -770 + Ag •77-k + A3 •77o7-k 



Direct comparison of (A3) and (A6) reveals 



Ai = A4 = 
A2 = -Ai - 
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Ae 
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-As , 
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(A6) 

(A7a) 
(A7b) 

(A7c) 

(A7d) 
(A7e) 
(A7f) 



These relations exhibit a redundancy in Bailin and Love's ansatz ( A4): only four of the eight 3- vector gap 
functions Ai, . . . , Ag are independent. The reason for this redundancy is that in Eq. (A4) the transverse 
components of A2, A4, Ag, and A7 never appear, as these gap functions are projected in the longitudinal 
direction k. Furthermore, only the transverse components of A3 and Ag appear on account of (Kq). Finally, 
the longitudinal components of the remaining two gap functions Ai and A5 can be absorbed by redefining 
A2 and Ag, thus only their transverse components contribute. In this manner, only half of the original 24 
gap functions are independent. Physically, this can be understood from the restriction to th e po sitive-parity 
channel J^ = 1+. Thus, only the four 3- vectors <fy^, TriJ: appear on the right-hand side of (A7). The other 
four 3-vectors 0_, 7r_ do not contribute, as they correspond to pairing in the J^ — 1^ channel. 
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This said, one can readily derive a more efficient form of Bailin and Love's ansatz (A4), which utihzes 
both longitudinal and transverse compone nts o f the t ruly independent gap functions. Choosing the latter to 
be Ai, A3, A5, and Ag, we obtain from (|a3) with (|A7|) 



A = Ai- (k + 7i(k)) +A3- (k + 7_L(k)) 7o7-k + A5- (k + 7i(k)) 70 + Ag • (k + 7i(k)) 7-k . 

(A8) 

It is now also easy to interpret the results obtained by Bailin and Love in ||2^, where they studied J^ = 1+ 
condensation in the ultrarelativistic limit. In this limit, a^ = 0, thus A5 = Ag = on account of (A_7). 
Bailin and Love discuss two order parameters for condensation in the J^ = 1+ channel, the first being the 
longitudinal component of A4 — A5 — Ag, the second being (Ai + A3) • (7 — k7 • k — ik x 775). 

From (A7), we identify the first with the longitudinal component of Ai + A3. From ( [A5| ), we realize 
that the second is identical to (Ai + A3) • 7j_(k) (1 — 707 • k), i.e., essentially the transverse components 
of Ai + A3. Thus, Bailin and Love discuss two separate gap equations, one for the longitudinal, the other 
for the transverse components of the independent gap function Ai + A3. Although we also find two gap 



equations for the longitudinal and the transverse gaps, cf. Section IV, we do not find that they decouple, as 
the excitation energies (127) contain both longitudinal as well as transverse gap functions. 
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